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Abstract 



We present new exact inhomogeneous vacuum cosmological solutions of Einstein's equations. They 
provide new information about the nature of general cosmological solutions to Einstein's equations. 
PACS 04.20.Jb, 98.80H, 04.30.-w 

Inhomogeneous and anisotropic solutions of Einstein's equations are of great interest because of the 
subtleties of their non-linearity and the clues they offer about the generic behaviour of cosmological models 
near spacetime singularities and at late times. They provide testing grounds for high-energy physics in 
the early universe and guide the development of quantum cosmology. Here, we present a new class of 
inhomogeneous vacuum solutions of Einstein's equations which generalise the Einstein- Rosen solutions [1]. 
They represent inhomogeneous counterparts of the most general Bianchi-type vacuum universes and may 
form a leading-order approximation to part of a general solution of the vacuum Einstein equations. Moreover, 
anisotropic universes containing perfect fluids with pressure less than the energy density behave in general 
like vacuum universes at early times. 

Consider the generalized Einstein-Rosen metrics: 

ds 2 = e^' z) (-df + dz 2 ) + g ab {i, z)dx a dx\ (1) 

where £ is a timelike coordinate; x a , a = 1,2, and z are spacelike coordinates. These metrics admit an abelian 
group of isometries, Gi- Spatially-homogeneous cosmologies of Bianchi types I- VII and the axisymmetric 
cases of Bianchi types VIII and IX, admit such a G2 and are particular cases [2]. Properties of (1) depend on 
whether = d^(det g a b)^ is spacelike, timelike or null (Greek indices run — > 3). The cases with a globally 
null or spacelike correspond to plane or cylindrical gravitational waves, respectively [2]. Metrics where 
the sign of B^B^ varies throughout the spacetime describe colliding gravitational waves [3] or cosmologies 
with timelike and spacelike singularities [4]. Metrics with a globally timelike B^ describe cosmological 
models with spacelike singularities. If the spacelike hypersurfaces are compact then the allowed spatial 
topologies [5] are a 3-torus, S 1 (g> S 1 (g> S 1 , for (detg a b)^ = £; a hypertorus, S 1 (g) S 2 , or a 3-sphere, 5 3 , for 
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{detgaiy)^ = sin z sin £ with < z, £ < ir. We shall present solutions for the globally timelike case. Metrics 
of this form have been studied approximately by Belinskii and Khalatnikov [6] who argued that they can 
provide the leading approximation to a general solution of the vacuum Einstein equations near the initial 
singularity. 

The metric g a b in (1) can be written as [6] 



9ab = R 



e a cosh 



sinh [3 



sinh (3 e 01 cosh (3 I 

where a = a(£, z), [3 = /?(£, z), and R = -R(£, z), so = R^. The vacuum Einstein equations are 

R-R" = 0, 



(2) 



(3) 



R R 1 
a + —a-a"- —a' = 2(a'/3' - a/3) tanh j3, 
R R 



22' 



+ \0 2 ~ P' 2 ) + \{# ~ «' 2 ) cosh 2 /3 + M (| ) - f | 



0, 



(4) 
(5) 

(6) 



R^_R 
R R 



(7) 



where dot denotes J| and a prime denotes J|. Physically, these equations describe the propagation of two 
inhomogeneous gravitational waves (a and /3) on an inhomogeneously anisotropic background spacetime. 
Assume that 



(3 = /(«), 

where / is to be determined. Eqns. (4) and (5) yield 



(8) 



d 2 f 



/2 ■ 2 

q — a 



= 0. 



(9) 



If a' 2 — a / 0, the general solution is 



2 



(3 = f(a) = artanh(^e a + Be~ a ) (10) 

with A and B constants. Using (9) in (4)-(5), and denning a new function F(u>) with u = u(£,z) by 
F(lo) = e a , we find two differential equations equivalent to (4) & (5): 

u> + & - J' - %rf = 0, (11) 
R ti 

d 2 F (dF\ 2 3 A 2 F 2 - B 2 F~ 2 - {1-2 AB) 

dJ + \dTj) (1 - 2AB)F - A 2 F S - B 2 F~ 1 ' ( ' 

Hence, for A / 0, 

e a = F{u) = [1 - 2AB - y/1 - 4AB tanhf^Vl - 4AB(a; - uj )}] I , (13) 
\/2A 

where ujq and fi > are constants, 4AB < 1, and tanh can be replaced by coth. 
For A = 0, F(uj) is found to be 

e a = F±(lo) = [B 2 + e ±2Mw-wo)]3. (14) 

If B = in (14) then /? = and the solution reduces to a diagonal metric of Einstein- Rosen form [1]. 

The behaviour of R(£, z) is needed to complete the solution and is determined by (3). We make a globally 
timelike choice for R to obtain cosmological solutions. Using the coordinate freedom, the general form is 

R{£,z)=Z. (15) 

With this choice the constraint equations for ip{S,,z) admit a first integral, and (6)-(7) become 

4 = -^ + ^[$ 2 + P' 2 + (a 2 + a' 2 ) cosh 2 /?], (16) 
$ = £\pfi + aa' cosh 2 0\. (17) 
Substituting F(u) from (13) or (14) into (16)-(17) yields 

ip' = W 2 &J, (19) 
where W 2 = fi 2 (l — 4AB), while equation (11) becomes 

Co + -u - oj" = 0. (20) 
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This is solved by 

oo 

w{£,z) =LO + a \ni + Tz + ^ cos[n(z - z n )]L4 n J «) + B n N Q (n£)], (21) 

n=l 

where A n , B n , loq, ao, z n , T are arbitrary constants (so we may set fx = 1); Jo and No are zero-order Bessel 
and Neumann functions. 

Define a function, G(£,z), which solves the following equations: 

G = ^ 2 W 2 ), (22) 
G' = (23) 
Eqns. (18) and (19) can be combined into a total differential 

dnj) = --d\n£ + W 2 dG 



and so 



z) = -hn£ + W 2 G(^ z) + M, (24) 



where M is a constant. 

Charach et al [7] have studied a system of equations very similar to (16)-(20) in a situation that corre- 
sponds to T = 0. Generalising to V / 0, we find from (22)-(23), 

G{t-,z) = lo + Ta z + ^a 2 In ^ + ^T 2 f 



+a cos[n(z - z n )][A n J (n£) + B n N (n£)} 

n=l 

oo 



n=l 

oo 



+W E n 2 ([A n J (nO + B n N Q {<)? + K-MO + B^nt)?) 



4 n=l 
^ oo 



ncos 2 [n(z - z n )]{A 2 n J (nOJi«) 



n=l 



+A n B n [N (nOJi(nO + J «)^iK)] + ^iVoK)iVi«)} 

-,00 oo 
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+ cos[n(z — z n )] cos[m(z — z. 



■in 



)}[mUM(0 - nff«(0]} 



where 



U^LiO = A n A m .h(ni)UnO + B n fi m JV K)iViK) + 2A n B m J 1 {ni)N {mi), 



(25) 



U^{0 = A n A m J {ni)J {mi) + B n B m N (n£)Ni (m£) + 24,B m J «)#i(m£). 



The full solution for the four unknown metric functions a, f3, ^,and R is given by eqns. (10), (13) or (14), (15), 
(21), (24), and (25). In order to understand the evolutionary behaviour of the solution we examine small and 
large values of n£. The initial curvature singularity of the universe lies at £ = and its expansion continues 
forever as £ — > oo. The product n£ can be interpreted as the ratio of the horizon distance in z-direction 
(Az = j d£) to some coordinate wavelength proportional to n _1 . The gravitational- wave inhomogeneities 
lie outside the horizon and evolve quasi-homogeneously when n£ < 1, but are eventually encompassed by 
the horizon and oscillate when n£ > 1, (density inhomogeneities behave similarly in non- vacuum models 
when the sound speed equals the speed of light, [8]). Here, we are interested in the overall behaviour of the 
cosmological model rather than that of particular modes of inhomogeneity and so n is assumed to be of 
order 1. We examine the solution as £ — > and £ — ► oo. It is convenient to express g a t in terms of F(u), 
since 




(26) 



The Limit n£ < 1 

There are two cases: A ^ and A = 0. 

(I) T/ie Case ,4 / : 

The asymptotic form of z) is 



oc 



cos[n(z - ,z n )]L4 n + -.B n ln n] 



n=l 

r> OO 



+ ln£[a + - Y] B n cos[n(z - z n )}\. 

IT A * 



n=l 



This implies that to — > — oo as £ — > 0. Hence, in this limit we have 



-»• -t^[1 - 2AB + Vl - 4AB]5 = F , 
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J AFo + BFq 1 

9ab U AFq + BF -i F -i 



Since 



^ 2 

G(£, z) — ► u; + Tao^ + oo > cos[n(z - z n )] [A n + -S n In n] 

i 7r 

n=l 

1 2 00 

+ ln^o«o + a o- 5 n cos[n(z - z n )]] 

2 7r — ' 

n=l 

we have -► e^Wf 012 ^ where 

oo 2 
cti(z) = M + W 2 ^o + W 2 Ta z + a W 2 V cos[n(z - z n )]L4 n + -S n ln n] (27) 

i 7r 

n=l 

W 2 1 2 °° 

"2(2) = ^r«o - o + «o^ 2 - 51 5 « cos i n ( z ~ z n)} (28) 
2 2 7r ' — ; 

n=l 

The change to proper time t is given by dt = e^d^ and so t ~ £ 2 . It is convenient to choose 0:2(2) 
such that the zero points of £ and t coincide; that is, with 02(2) > —2. We see that g a b has off-diagonal 
components and is linear in £ as £ — ► 0. As £ — ► the universe shrinks along the x- and y-axes, but the 
behaviour along the z-axis depends on the overall sign of 02(2). In the range — 2 < 02(2) < the expansion 
along the z-axis blows up and produces a spindle singularity. However, choosing 02(2) to be positive results 
in a universe shrinking in all direction to a point-like singularity In general, the oscillatory character of 02(2) 
will lead to a complicated oscillatory spatial structure with some regions expanding along the z-direction 
whilst others implode according to the relative magnitudes of the constants uq, W, and B n . 
(II) The Case A = : 

The interesting case arises from F = F- in (14) since F + — ► B as uj — > — 00, which reduces to case (I). 
However, F_ — > 00 as u; — ► — 00 with F_(lj) ~ £- a s(*) where 

2 00 



o 3 (z) = ^[oo + ~^2 B n cos[n(z - Zn)]] (29) 

and so 



n=l 



9ab ~ I B ^ 1+a3 ( z ) £l+a 3 (z) 
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Neglecting the off-diagonal components as £ — > 0, and re-expressing £ in terms of proper time t, in regions 
where 0:3(2) > the approach to the singularity has the form 

2(1-0:3) 2(1+0,3) 2ct 2 (z) 
(t 2Pl ,t 2P2 ,t 2p3 ) = (t«2(«)+a,ta2M+2,t«»2M+2), 



so EtiP< = i and E*=iP? > i- 



The Limit n£ > 1 
(I) The Case A ^ : 

The asymptotic form of a> is 



— > ujq + Tz + q In £ 

S ( — ) cos[n(z - * n )] [A n cos« - y ) + B n sin« - 7)] (30) 
_ \ Tin J 4 4 



^ \im) 

so — > 00 as £ — > 00. Now 



so that 



G(£,z) -> u; + ro z+ -Oohi£ 

n=l 



J F AFq + BFq 1 

9ab^U AFo + BF -i F -i 



and — > e 71 ^^ 72 e 73 ^ +74 ^ 2 as £ — > 00, where the four 7, are defined by 



{71(2), 72} = {M + W 2 Wo + ^ 2 ra z, ^W 2 a§ - i}, (31) 



00 W 2 
{73, 74} = {^E ^ + B nh — ^} (32) 
n=l 
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The £(i) relation cannot generally be integrated in this case since 



72 
2 



e 2 d£ 



However, if the choice 72 = 2 is made then t ~ exp{274(£ + A) } + to, to constant, where A = 73/(274) . 
If 74 > 0, t is a monotonically growing function of £ for £ > 0, so in this particular case the behaviour of 
the solution in £ time reflects the behaviour in proper time t. If T = 74 = then t ~ £T 2 / 2 exp{73^/2}. In 
the limit £ — > 00 anisotropic expansion is found: the expansion rates along the x- and y-axes equalise but 
differ from that in the z-direction. There can be no genuine isotropic limit because there is no flat vacuum 
Friedmann universe. 

(II) The Case A = : 

We have F_(lj) — > B as to — > 00 (this is equivalent to Case (I)), whereas F+(u) — > 00. Neglecting the 
oscillatory part of w(£, z), the asymptotic behaviour of F + is -F+(o;) ~ £^ a °, and hence 

(£l+/nao ^jtl-juao \ 

Recalling that ~ £ 72 exp{73^ + 74^}, and neglecting the off-diagonal components of g a b, we see that this 
describes an anisotropic universe at late times. 

Our solutions are of particular interest for the problem of determining the meaning and characterisation 
of a general solution of the Einstein equations near a cosmological singularity [9] . It was argued in [6] that a 
metric of the form (l)-(2) is a leading-order approximation to a 'general' cosmological solution of Einstein's 
equations characterised by four independent arbitrary functions of three space variables on a spacelike hy- 
persurface of constant time. This requires that derivatives with respect to the variables x a to be smaller 
than those in the z and t variables, so that the x a appear only as functions of integration, and that the 
metric component ^(omitted from the metric (l)-(2)) be smaller than g zz . These requirements are met if 
9zz/9ab *C 1, [6]. When this condition holds our solution lays claim to being a leading-order approximation 
to part of a general cosmological solution of Einstein's equations [10]. It is an inhomogeneous generalisation 
of a single cycle of Mixmaster oscillations during which one axis (g zz ) changes monotonically whilst the other 
two oscillate [11]. Now, in the limit n£ <C 1, when A ^ we have g zz /g ab ~ f* 2 " 1 ~ t 2 ("2-i)/(a 2 +2) ^ 1 
only if a 2 > 1; however, when A = we have g zz /g xx ~ £"2+^3-1 an d g zz /g xy ~ g zz /g yy ~ g<*2-a 3 -i^ so 
g zz /g a b *C 1 requires 02 > 1 and [03) < a.2 — 1. The on are defined by (27)-(29). When n( > 1 we have, 
for A / 0, that g zz /g a b ~ £ 72-1 exp{73£ + 74^ 2 }. The 7$ are defined by (31)-(32). Since 73 and 74 are 
positive we have g zz /g a b 3> 1 in this case, and our solution is not a leading-order approximation to the 
general solution when £ — > 00. When A = for n£ » 1, we have g zz /g xx ~ £T2-i-Mao exp{73^ + 74^} and 
gzz/gxy ~ £T2-i+Mao e xp{73^ + 74^ 2 }- Again, as in the ^4/0 case, the solution cannot be the leading-order 
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approximation to a general (4-function) solution. Thus for certain parameter choices and over certain in- 
tervals of time our solutions provide new information about a general solution of Einstein's equations near 
a cosmological singularity. 
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